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Abstract 

We study electroweak baryogenesis and electric dipole moments in the presence of the two 
leading-order, non-renormalizable operators in the Higgs sector of the MSSM. Significant qual- 
itative and quantitative differences from MSSM baryogenesis arise due to the presence of new 
CP-violating phases and to the relaxation of constraints on the super symmetric spectrum (in par- 
ticular, both stops can be light). We find: (1) spontaneous baryogenesis, driven by a change in 
the phase of the Higgs vevs across the bubble wall, becomes possible; (2) the top and stop CP- 
violating sources can become effective; (3) baryogenesis is viable in larger parts of parameter space, 
alleviating the well-known fine-tuning associated with MSSM baryogenesis. Nevertheless, electric 
dipole moments should be measured if experimental sensitivities are improved by about one order 
of magnitude. 
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I. INTRODUCTION 



Electroweak baryogenesis (EWBG) is an attractive mechanism for generating the baryon 
asymmetry of the Universe (BAU). Its primary attraction is the possibility to experimentally 
test two of the three Sakharov conditions. Aspects of the departure from thermal equilibrium 
via the electroweak phase transition (EWPT) can be explored in collider experiments, while 
CP violation can be tested in electric dipole moment (EDM) searches. 

In the EWBG picture, electroweak symmetry breaking proceeds via a first-order phase 
transition, where bubbles of broken SU(2)x, symmetry nucleate and expand in a background 
of unbroken symmetry. CP-violating interactions within the bubble walls lead to the pro- 
duction of CP-asymmetric charge density of left-handed fermions. This charge, diffusing 
ahead of the wall into the unbroken phase, is converted into the BAU by non-perturbative, 
electroweak sphaleron processes. To the extent that electroweak sphalerons are inactive after 
electroweak symmetry breaking, the baryon density "freezes out" once it is captured by the 
advancing bubble wall. This mechanism satisfies the Sakharov criteria [l| and generates the 
BAU provided two conditions are met: (1) the phase transition is "strongly" first-order (oth- 
erwise electroweak sphalerons are active within the broken phase and washout the BAU), and 
(2) the CP violation is sufficient to generate the observed BAU. Neither of these conditions 
are met in the Standard Model (SM) [^|. 

Beyond the SM, the most widely studied EWBG model has been the Minimal Super- 
symmetric Standard Model (MSSM). However, there are a series of tensions that make 
this scenario severely constrained by experiment. First, there is tension in the top squark 
("stop") sector. A strong first-order phase transition requires at least one light stop (which 
must be mostly tn, to avoid large contributions to the p parameter and due to null searches 
for a light sbottom js, 4]). At the same time, the large radiative corrections needed to push 
the Higgs boson mass above the LEP bound > 114 GeV (4] require that at least one 
stop (ti) is very heavy Recently, the phase transition was studied in an effective theory 
with a large stop hierarchy, concluding that successful EWBG is possible only for mj R < 125 
GeV and > 6.5 TeV |6J, rendering the scenario finely tuned. 

Second, there is tension between having enough CP violation to produce the BAU and 
evading stringent constraints from EDM searches. In the MSSM, the CP-violating phases 
that drive EWBG arise in the gaugino/higgsino sector. The same phases contribute to 
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EDMs. While one-loop contributions can be sufficiently suppressed by making the first two 
squark and slepton generations heavy, there exist two-loop contributions that cannot be 
suppressed without spoiling EWBG (assuming no fine-tuned cancellation between different 
EDM contributions) and which predict a minimum value of the EDM. These "irreducible" 
EDMs strongly constrain the viable MSSM parameter space: EWBG with universal gaugino 
phases is nearly ruled out. With improvements by a factor 3 — 4 in the upper bounds on 
the EDMs of the electron or the neutron, MSSM baryogenesis will be possible only in the 
so-called "bino-driven" scenario, where the CP-violating phase associated with the U(l)y 
gaugino is tuned to be much larger than that of the SU(2) L gaugino [7, 8]. 

Third, there is tension in the mass of the pseudoscalar Higgs boson A . Large values of ijia 
are preferred (i) to make the EWPT more strongly first-order, and (ii) to evade constraints 
from b — > S7 [8|. However, the production of left-handed charge during EWBG is enhanced 
when rriA is light. There is also a tension in the value of tan/3 (i) from a compromise in 
giving a large enough value of the Higgs mass versus a strong enough phase transition, and 
(ii) from the constraints from b — > for small values of tua- All in all, from the theoretical 
point of view, these tensions force the MSSM (if it is to account for EWBG) into a narrow, 
finely tuned region of parameter space. 

An attractive extension of the MSSM is the "Beyond the MSSM" (BMSSM) scenario Q. 
Here, a non-renormalizable contribution to the MSSM superpotential is included 

H^BMSSM = H-MSSM + "TJ (H u H d ) 2 , (1) 

as well as a contribution to the soft SUSY-breaking Lagrangian 

/-BMSSM /-MSSM , A s ^SUSY / rr tt \2 f n \ 
L soft = ^soft H u d > ' > 

encoding the leading supersymmetric and F-term supersymmetry breaking corrections to 



the Higgs sector that arise from a new threshold at mass scale M 10Ml4l] . The corrections 



enter the spectrum and interactions through the dimensionless parameters 



For M ~ few TeV, the BMSSM has interesting implications for cosmology 
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and for 



Higgs phenomenology [20J. The BMSSM operators, which contribute at tree-level to the 
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TABLE I: The CP-violating phases in the BMSSM. Here i = 1, 2, 3 labels gaugino mass parameters 
and f labels the trilinear sfermion-Higgs coupling corresponding to a SM fermion f. 

Higgs mass, alleviate the tension associated with the stop sector and tan (5. Now, the left- 
handed stop can also be relatively light, providing additional bosonic degrees of freedom 
that strengthen the first-order phase transition 1 . 

In this work, we examine the BMSSM implications for CP violation and the generation of 
the baryon asymmetry. In Sec. [Til we describe new CP-violating phases associated with the 
BMSSM operators. In Sec. EEl we review relevant aspects of the phase transition dynamics 
and show that these phases lead to new CP-violating sources that generate charge density 
during the phase transition. In this section we also compute the resulting BAU. In Sec. IIVI 
we discuss how searches for EDMs constrain CP violation and baryogenesis in the BMSSM. 
We conclude in Sec. [V] The appendices contain details of the CP- violating vacuum structure 
and radiative corrections to the Higgs CP-violating phase. 



II. CP VIOLATION 

In this section, we describe the BMSSM Lagrangian to leading order in M _1 , emphasizing 
those aspects that are relevant for CP violation, baryogenesis, and EDMs. The new BMSSM 
phases (denoted $1,2) lead to (i) explicit CP violation in the neutralino, chargino, and squark 
mass matrices, and (ii) CP- violating mixing of the Higgs pseudoscalar Aq with the two other 
neutral Higgs scalars ho,H . We express our results in terms of physical CP- violating phases 
that are invariant under phase redefinitions of fields, summarized in Table [U 

Neglecting flavor mixing, the invariant phases of Table [TJ provide a complete basis for all 
of the rephasing invariants in the BMSSM. Such basis is easily constructed by noting that 



Higgs field rephasing is equivalent to global U(1)pq and U(1)r_pq transformations 22 ] . 



1 To be clear, a first-order phase transition is induced radiatively through thermal effects, similar to the 



MSSM, as opposed to new tree-level interactions as in, e.g., the NMSSM [21 1. 
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The U(1)pq and U(1)r_pq are explicitly broken by the dimensionful MSSM parameters 
appearing in Table HI as well as by the BMSSM new effective couplings. By promoting 
the parameters to spurions with well-defined transformation properties, one can extract the 
rephasing invariants in terms of U[1)pq and U(1)r_pq conserving combinations. 



We follow the notation of Ref. 



23[ with respect to the MSSM parameters. In our numer- 



ical analysis we implement the quantum corrections from the neutralino, chargino, scalar 
Higgs, and squark sectors. Details are given in Appendix [Bj CP violation induced by these 
corrections in the Higgs sector is suppressed for small values of the trilinear A term and for 
moderate values of /i, which we adopt throughout our analysis. This allows us to focus on 
the novel tree-level BMSSM effects. 

First, we consider the tree- level Higgs potential 

, 2 , g' 2 + 9 2 



V = (m 2 Hu + |/i| 2 ) \H U \ 2 + [m\ d + \n\ 2 ) \H d 



2\ 2 



\H U \ — \H d \ ) 



9_ 
2 



+ (b{H u H d ) + 2e l (\H u \ 2 + \H d \ 2 ) (H u H d ) + e 2 {H u H d f + h.c. 



(4) 



with SU(2)x, contractions defined as (H u H d ) = H+H d — H^H® . At zero temperature, the 
Higgs vacuum expectation values (vevs) are 



'X) 



sp v e 



\ H d) = v d = c p ve 



(5) 



where tan/3 = \v u /v d \, s@ = sin/3, cp = cos/3, and v ~ 174 GeV. The relative phase (9 u —9 d ) 
is unphysical and can be set to zero by a gauge transformation. We define the Higgs phase 
9 as 



e u + 9 d = 9- arg(6) . 



(6) 



It is useful to factor out arg(6) explicitly, since 9 is rephasing invariant 24J 2 . Next, we 
define rephasing invariant BMSSM parameters 



€i r = \ei\ cos($i + 9) 
e 2r = \e 2 \ cos($ 2 + 26*) 



en = |ei| sin(^i + 9) 
e 2i = |e 2 | sin(^ 2 + 29) 



(7a) 
(7b) 



2 In Ref. [25|, explicit CP violation in the BMSSM was also studied. However, the authors considered a 
scenario in which arg(6) = 8 U + 6d = 0. Since 9 is rephasing invariant, Ref. [25j deals with a very specific 
physical model, and the results derived there do not apply in general. 



5 



In the phase convention 6 U + B& = 0, our definitions reduce to the usual definitions e\ r = 
Re[ei], etc. (a]. 

The masses and mixing angles of the physical Higgs bosons receive tree-level corrections 
from the BMSSM operators. In our expressions to follow, we work at tree-level and eliminate 
the set of parameters (m 2 Hu ,m 2 ld , \b\) in favor of (v , tan/3, m^). We parametrize the Higgs 
fields in the following way: 



H,„ 



Q r> 1 1 J- ^ iu ~^~'^ au 



Ha 



(8) 



In the limit en = 
as in the MSSM 



e 2 



0, one can separately diagonalize the CP-even and odd Higgs states, 



231 ] . The eigenstates are 







Uo j 





cos a — sin a 
sin a cos a 





sin /3 — cos /3 ' 


II 




cos /3 sin /3 J 




[a d J 



(9) 



with Higgs mixing angle a given by 



cos 2a 



m| + 4e 2r w 2 



m 



_ m h 

The mass eigenvalues also receive tree-level contributions proportional to e% t 2r 
ticular, the lightest Higgs boson receives a correction 



cos 2/3 



sin 2a 



{m 2 A + m 2 z ) sin 2(3 — 8ei r w 2 



m 



mi 



■ (10) 



9|. In par- 



8 t m 2 h 



2e lr {m\ + m 2 z )s2f} + e 2r (mi - m 2 z )c 2 w 

2v I e 2r - 2 ei r s 2/3 

2 „2 



— m|) 2 + 4m 2 A m 2 z s 



2rf 



This contribution can increase the tree-level Higgs mass above the LEP bound, without the 



need for radiative corrections 26]. It is important to recall that the LEP bound on the 
ightest neutral Higgs boson mass is drastically changed in the presence of CP violation 
27| . In fact, there are allowed regions even for very small values of the Higgs boson mass. 
In principle this implies that a much larger region in parameter space can accommodate a 
strong first order phase transition and the right-handed stop can be heavier than the top 
quark. However, as we will see below, there are significant constraints on the amount of CP 
violation. 

CP violation enters the Higgs sector at tree-level when e\2% 7^ 0, leading to mixing between 
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CP eigenstates. In the (ho, Hq, Aq) basis, the Higgs mass matrix is 



m 2 h m\ A ^ 
m 2 H m 2 HA 



(12) 



\ m hA m HA m A ) 

The remaining CP-odd state Go is eaten by the Z boson. The parameters m\ A and m 2 HA , 
which govern the mixing between CP-even and odd states, are given by 

m\ A = 4v 2 eu sin(/3 - a) - 2v 2 e 2i cos(a + /3) « -2v 2 (e 2 » s 2/3 - 2e Xi ) , (13) 
m HA = 4w 2 eii cos(/3 - a) - 2v 2 e 2i sin(a + /?) m 2v 2 e 2i c 2j3 , (14) 

where the approximations follow in the limit of moderate tan/3 and m\ 3> m 2 z , such that 
a w /3 — 7r/2. To 0(ei,2i), the eigenvalues are unchanged from the CP-conserving case. Note 
that ?Ti is now the mass of the "mostly-pseudoscalar" eigenstate, not the mass of Aq. To 
avoid this confusion, we will express physical quantities in terms of the charged Higgs boson 
mass, using the relation 

m 2 H± = m 2 A + m 2 w + 2e 2r v 2 . (15) 

The Higgs phase 6, determined by the minimization condition dVo/dO = 0, is given by 

, a 2v 2 (e 2i s 2fS - 2 eu) , . 

tan^ = — r . (16) 

S2fiK ± -™>w) + 2v 2 {e 2r s 2/3 - 2e lr ) 

In the small tan (3 regime (such that cot f3 ^> \e\\v 2 / m 2 A ) , one can treat 9 perturbatively since 
it is 0(ei j2 ) ■ However, in the large tan/3 regime (such that cot /3 < \e\\v 2 /m 2 A ) , one can 
have = 0(1). In this regime, the Higgs potential can develop more than one minimum 
in the direction; we discuss this possibility in Appendix [A] In practice, constraints on 



b — > s 7 (28| imply that the mass of the charged Higgs cannot be too light (rriH± ^ 300 GeV) 
unless the charged Higgs contribution to b — > 57 interferes destructively with some other 
process 3 . Restricting ourselves to > 200 GeV and tan/3 < 10 is sufficient to avoid 

additional phase minima. 



3 Some amount of interference is in fact expected, considering the light stops and charginos of our baryo- 
genesis scenario, weakening the bound on tuh± ■ 
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(19) 



CP violation from the complex BMSSM parameters ei )2 also enters the SUSY mass ma- 
trices, potentially impacting both baryogenesis and EDMs. For example, the top squark 
mass matrix, in the (t* L , t* R ) basis, is 

m 2 = I m Q 3 + V? S W + yt(A*v* u - fiv d + 2e lVu v 2 d /fi*) \ ^ 

with D-term contributions = (T^ — sin 2 9w) C2/3 m 2 z . We define the stop mixing 
parameter X t = | [m~] 12 | as the magnitude of the off-diagonal entry in Eq. (1171) . The chargino 
mass matrix, in the (W + , H+, W~ , H d ) basis, is 

/ X T \ I M 2 gv* u 

M 5 =[ , X = , I ■ ( 18 ) 

yX J \gv* d fi - 2e 1 v d v u / fi* 

The neutralino mass matrix, in the (B, W°, H d , H®) basis, is 

( M x -9'v* d /y/2 g'v*jV2 \ 

M 2 gv* d /V2 -gv*/V2 

-g'v* d /V2 gv* d /V2 2e x v 2 Jii f -fi + Aeiv u v d /fj,* 

V g'v* u /V2 -gv* u /V2 -12 + Ae x v u v d /n* 2e 1 v 2 / j j* ) 

In each of these mass matrices, the BMSSM parameters lead to new sources of CP violation 
through both the explicit factors of e x and the complex phase of the Higgs vevs, depending 
on ei and e 2 . Physical CP violation observables depend on these phases only through the 
invariant combinations listed in Table HI 

To summarize, CP violation from the BMSSM manifests in the following ways: 

• Mixing arises between CP-even and CP-odd neutral Higgs eigenstates, proportional 
to rn} hA and m 2 HA . 

• The Higgs phase 9 enters through the Higgs vevs in the SUSY mass matrices. 

• The parameter e\ appears explicitly in the SUSY mass matrices. 

In Sec. PH we show that the BAU induced through BMSSM phases is approximately pro- 
portional to 9. In Sec. IIVI we find that the dominant contributions to EDMs arise through 
either 9 or Hq-Aq mixing. Since 9 and vn? hA are proportional to the same linear combination of 
61,21, EDM constraints will provide direct bounds on the EWBG mechanism in the BMSSM. 
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III. ELECTROWEAK BARYOGENESIS 



In this section, we describe how electroweak baryogenesis is realized in the BMSSM. First, 
we study the nature of the phase transition and the properties of the expanding bubbles 
relevant for the BAU computation. Second, we identify novel sources of CP violation in 
the BMSSM and compute the resulting BAU. The new sources are induced by a variation 
of the Higgs phase 9 across the bubble wall. Therefore, we devote special attention to the 
computation of the temperature and space-time dependence of 9. Similar effects arise also 



in the MSSM at the quantum level 29, 



30|. In the BMSSM they arise classically [l5|, and 



can be quantitatively much more significant. 



A. Phase Transition and Bubble Properties 

1. The critical vev and temperature 

The EWBG mechanism requires a "strong" first-order phase transition to avoid sphaleron 



erasure of the BAU within the broken phase. This condition is satisfied if 31] 

^>l, (20) 

where T c is the critical temperature (defined here as the temperature of free-energy degen- 
eracy between the broken and symmetric phases) and v c = v (T c ) is the Higgs vev in the 
broken phase at T c , in the normalization of Eq. ([S]). Since a first-order phase transition in 
the BMSSM arises through radiative corrections involving stops (as in the MSSM j^, 32]). 



Eq. (1201 provides important constraints on the parameters of the stop sector. In addition, 
the size of v c itself is important for the BAU computation; as we show below, the BAU scales 



as v\. 



We compute v c and T c using the two-loop finite-temperature effective potential of Ref. |11 
provisionally neglecting the effect of CP violation. In Fig. [JJ (left panel), we show how T c 
and v c depend on the stop parameters [c.f. Eq. (fl7|) ]. We consider two cases: mixing 
(X t = (150 GeV) 2 ) and no mixing (X t = 0), while varying mfj 3 (assuming m^ 3 < 0) and 
fixing rriA = 250 GeV, tan/3 = 5, = 114 GeV, and mq z = 200 GeV. At a given value of 
mfj 3 , stop mixing suppresses v c /T c , thereby weakening the phase transition. The filled circle 
at the lower edge of the X t > line (gray) corresponds to the maximal value of m^ 3 where 
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FIG. 1: Left: The critical vev v c (solid curve) and temperature T c (dashed) vs. the right-handed 
stop soft mass mjj 3 (imaginary, since mfj 3 < 0), for two values of the mixing parameter Xt ■ Right, 
black: The minimum value of mg 3 , consistent with the electroweak oblique parameter T < 0.2, as 
a function of the stop mixing \fX~t, for two values o/tan/3. Gray: Model independent, direct limit 
from sbottom search. 

Eq. ( 120]) is fulfilled. On the other hand, decreasing m^ 3 increases v c /T c , strengthening the 
phase transition; eventually, however, this leads to an undesirable tachyonic stop, denoted 
by the filled circle at the upper edge of the X t = line (black). These effects are well 
known in the MSSM: stop mixing effectively screens the one-loop cubic correction to the 
effective potential and increases the value of the Higgs mass, while such screening can be 
compensated by a sufficiently negative m^. 

Next, we consider how the strength of the phase transition depends on rriQ 3 . We vary rriQ 3 
while keeping the neutral Higgs mass fixed by simultaneously adjusting e li2r . For rriQ 3 > 200 
GeV, v c and T c are only weakly dependent on mg 3 . For vtlq 3 < 200 GeV, the phase transition 
is strengthened as one decreases m,Q 3 , allowing for Eq. (J20l) to be consistent with greater 



values of mjj 3 



lol . \r\ . However, experimental constraints provide a lower bound on mg 3 , 
shown in Fig. [1] (right panel). The collider search bound on bottom squarks (m^ > 89 
4|) implies that mg 3 > 80 GeV 17|, shown in gray 4 . The precision constraint on 



GeV 



the p parameter, which is sensitive to the 0(rn^/mQ 3 ) relative mass splitting in the doublet 



For a neutralino LSP with mass m x o < 90 GeV, the bound is significantly stronger: > 250 GeV 33 1. 
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{pLi provides a stronger bound [3|,|4j. We find that rriQ a > 160 GeV is required at the 95% 
CL, unless the stops are substantially mixed; at the 90% CL one must keep mg 3 > 200 GeV. 



2. The bubble profiles 



We now turn to the properties of the bubbles that nucleate and expand during the phase 



transition 



34 



35]. These bubbles are characterized by spacetime-dependent background 
Higgs fields, described by (v, (3, 9). In the MSSM, the dominant sources driving EWBG are 
proportional to the small parameter A/3 < 10 -2 associated with the spacetime variation of 
(3 34j. Here, we neglect variation in /3 (A/3 = 0) in order to focus on the leading BMSSM 
effects proportional to the analogous parameter A9 associated with variation of 9 across the 
bubble wall. 

The bubble profiles for v and 9 are determined by equations of motion 



dV T 
dv 
dV T 



d 2 v(z) 



dz 2 



'2/3 



89 2 

assuming a vanishing Z background 30, 



d z (v 2 (z)d z 9(z)) , 



(21a) 
(21b) 



Here, Vr is the finite temperature effective 



potential and the coordinate z is the distance from the wall in its rest frame (we assume a 
planar bubble). The boundary conditions are such that z = +oo (— oo) corresponds to the 
broken (unbroken) minimum of the potential. 

Rather than solving Eq. ( 12 lap , we assume a kink ansatz for the Higgs vev: 



v{z) 



v(T) 



tanh 



2 z 



(22) 



To avoid complication, we set v(T) = v c . This potentially underestimates the BAU: since the 
time of bubble growth necessarily corresponds to T < T c , one in general expects v(T) > v c , 
resulting with more effective sources. The approximation is justified if, following the onset 



of the phase transition, the universe is reheated back near the critical temperature 



33, 



The wall width L w is defined to match the kink ansatz of Eq. (1221) onto the bubble profile, 

(23) 



where v m i n = 0.1 v c and t> max = 0.9 v c designate the field on either side of the wall. Neglecting 
CP violation in Eq. ( 123]) . we find values in the range L w = (15 — 40)/T c . 
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Next, we obtain the profile of the vev phase 0(z). In order to clearly illustrate the essential 
dynamics, we make two simplifications: (i) we consider the < 1 regime, satisfied when 
cot /3 ^> \ei\v 2 /m\, and (ii) we neglect all radiative corrections to the phase-dependent part 
of Vt- Under these assumptions, Eq. f!21bj) becomes 

v 2 (z) 



d z (v 2 (z)d z 9(z) 



v 0) (%, - %) 



6{z) 







(24) 



In the present section, the zero temperature Higgs phase and vev are denoted 9q and Vq for 
clarity Eq. ( 1241) can be cast in the dimensionless form 



dr(gd r f) = rg(f-g) 



where 



/(*) 



v 2 e 



2z 



[m 



T 



m 



Wi 



9{r) 



1 + tanh r 



(25) 



(26) 



For m# ± > 200 GeV, we have r > 100. Thus, the solution to Eq. (j25|) is governed by the 
potential energy term on the RH side, such that f(r) ~ g(r) + 0(t~ 1 ). This conclusion 
is borne out by numerical evaluations for r > 30, which we perform using the method of 
Ref. 



In summary, the Higgs phase profile is given by 



9{z) 



A9 



tanh 



2 z 



(27) 



where m 9qv 2 /v 2 . We note that 9(z) is proportional to the square of the kink in Eq. (I2~2l . 
not linear. Furthermore, since EDMs are directly sensitive to the value of 8q, this sensitivity 
translates into a direct constraint on the phase variation across the bubble wall, and hence 
on baryogenesis. 

The preceeding analysis can be generalized away from the 0^1 regime, necessary when 
tan/3 > 10. Assuming that the profile of 9 is again dominated by potential energy (such 
that dVr/d0 ~ 0), we find the following approximate solution for 9(z): 

2(s 2 /3 e 2i -2e lt ) {v(z) 2 - ug) 



tan (9(z) - 9 



'2/3 



(m 2 H - m 2 w ) + 4 e lr (v(z) 2 - v 2 ) ' 



(28) 



This solution gives the leading behavior of 9(z) in all tan (3 regimes: (i) it reduces to Eq. ()2j 
when cot/3 3> \ei\v 2 /m 2 A limit, and (ii) it is valid to leading order in cot/3 when cot/3 < 
\e\\v 2 /m 2 A . 
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In Appendix [B], we study the impact of radiative corrections on the Higgs phase at zero 
and finite temperature. In particular, we find that 6(z) is shifted by an overall constant, 
while A6 remains approximately unchanged. 



3. Wall velocity 



The bubble wall velocity is an important parameter in the EWBG computation. A recent 



0.4 in the MSSM 



401 ]. significantly larger than previous estimates of 



study found tha t v v 

v w ~ 0.01 — 0.1 41]. Therefore it is worthwhile examining how the BAU depends on v w . 

The optimal wall velocity for EWBG arises as a competition between two Sakharov 
conditions. The generation of baryon number (ub) is fueled by chiral charge diffusing ahead 
of the advancing bubble wall, characterized by an effective diffusion constant D and a 
diffusion time r^g = D/v^ 42J. If electroweak sphalerons are in equilibrium, with rate 



T W s 3> Tdiff' n B is suppressed, as per the third Sakharov condition. On the other hand, if 
T ws <C t^, then ub is also suppressed, since few baryon number violating processes occur. 
Therefore, the maximum baryon number production occurs when T ws ~ r^g, corresponding 
to a velocity v w ~ a/ DT ws ~ (few) x 10~ 2 38]. In our numerical computation, described 
below, we indeed find that ub is maximized for v w = 0.03. 

If we consider the range 0.01 < v w < 0.4, we find tib varies by a factor of 4 — 5, with the 
minimum ub for v w = 0.4. For the sake of definiteness, we fix v w — 0.1, which is approx- 
imately the central value for n#. We expect that v w in the BMSSM can be approximated 
by the MSSM case. Potentially, the presence of the light LH stop leads to an additional 
contribution to the frictional force determining v w . However, we expect this to be a minor 
effect since ti, which cannot be too light, is somewhat Boltzmann suppressed and acquires 
only a fraction of its mass via the Higgs mechanism. 



B. Baryon Asymmetry Computation 

The computation of the BAU involves a system of coupled Boltzmann equations of the 
form 

d t n a - D a V 2 n a = T ab n b + Sf . (29) 

b 
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Here n a is the charge density for species a. The CP-violating source Sfj^, non-zero only within 
the moving bubble wall, leads to the generation of non-zero n a . The diffusion constant D a 
describes how efficiently n a is transported ahead of the wall into the unbroken phase where 
sphalerons are active. The interaction coefficients r a ;, correspond to (i) inelastic processes 
that convert charge from one species to another, and (ii) relaxation processes that wash out 
charge within the broken phase. Although BMSSM contributions modify T a b at 0(ei^), it 
is safe to neglect these corrections. Previous studies have shown that the solutions to the 
Boltzmann equations are insensitive to sub- 0(1) variations in the interaction coefficients 



42 



441 ] . We refer the reader to Ref. 43|, |44j, which we follow here, for details concerning the 



setup and derivation of the Boltzmann equations in the MSSM. 



1. CP-violating sources 



The novelty of BMSSM baryogenesis appears in the CP-violating sources. We compute 



45 



46l | . More sophisticated 



47 



48]. 
49]. 



these sources following the "vev-insertion" approach of Refs. 
treatments, going beyond the vev-insertion approximation, exist in the literature 
However, there remains some controversy, and this is an area of active investigation 
Therefore, we opt for the simplest framework (vev-insertion) by which we may present our 
new BMSSM sources. 

The higgsino CP-violating source, which drives EWBG in the MSSM, receives important 
BMSSM contributions. We have 



3g 2 K 2 (z) 
2vr 2 



k 2 dk 



Im 
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+ 
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Im 



1 - n F (£g) - n F (£g) 



(30) 



The first and second terms correspond to the sources induced through higgsino-wino and 
higgsino-bino mixing, respectively. The important BMSSM effects enter into the prefactors 



Ki{z) = \M if i\ v\z) sin (fa + 9(z)) $(z) + ^ cos (0, + 9(z)) 9(z 



(31) 



where M 12 are the gaugino mass parameters. In the MSSM we have 9 = 9 = 0, so that the 
sources are driven by the gaugino phases fa defined in Table HI However, in the BMSSM, 
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contributions arise from ei^i, entering through 9. Futhermore, the second term in Eq. (1511) . 
which is unique to the BMSSM, is not suppressed by A/3 < 10 -2 . The momentum integrals 
in Eq. (I50|) are identical to the MSSM case, discussed in Ref. [46|; roughly speaking, they are 
maximized "on- resonance" (when |M,| ~ and are highly suppressed far off-resonance. 

In the MSSM, the CP-violating sources for third generation squarks cannot drive EWBG. 
The Higgs mass bound requires that and bi are heavy, and thereby Boltzmann suppressed 
in the electroweak plasma. The BMSSM opens the door for squark-driven baryogenesis, since 
the stops and sbottoms can be relatively light. The CP- violating sources for stops (q = t) 
and sbottoms (q = b) are 



s 



(32) 
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2tt 2 

with prefactor 

Kq(z) = \A qf i\ v 2 (z)P(z) sin(<p q + 6{z)) 

+ V ~jt (stp \A ql i\cos^ q + 6{z)) + s% (\fi\ 2 -\A t \ 2 )) 9(z). (33) 

In addition to the squark phases <p q appearing in Table [J, the CP- violating sources include 
contributions from ex s 2i that enter through 9. 

In the BMSSM, there are CP- violating sources for the third generation quarks, top (q = t) 
and bottom (q = b), due to the Higgs phase 9(z), that does not arise in the MSSM: 
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where 



K q {z) 



v 2 (z) s%9\z). 



(34) 



Again, we refer to Ref. 



46] for the notation of quantities within the momentum integral. 



We note that because the quark thermal masses and widths are approximately equal for 
qi and q^ (dominated by common QCD effects), the momentum integral is suppressed. 
This situation may be an artifact of the vev-insertion approach. Similar quark CP- violating 
sources, computed within the WKB approximation, have been studied within the contexts 



of Two Higgs Doublet models [5Cj and the SM with higher dimensional operators 



5l| 
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FIG. 2: The baryon asymmetry generated by Q -terms in CP-violating sources in units of the observed 
BAU. Left panel: top source (blue, dashed), stop source (green, dashed), and their sum (black, 
solid) as a function of the soft LH stop mass parameter. Right panel: higgsino-bino source (purple, 
dashed), higgsino-wino source (red, dashed), and their sum (black, solid) as a function of \i. Other 
relevant parameters are specified in Table O 

2. The baryon asymmetry 

Here, we compute the baryon asymmetry. As argued above, we expect that the leading 
BMSSM effects will enter through the CP-violating sources. We neglect BMSSM effects 
arising in the various transport coefficients and diffusion constants that enter the Boltz- 
mann equations, following the general MSSM setup described in Ref. 44J. In addition, we 
make the further assumption of chemical equilibrium between particles and their superpart- 
ners, valid when gauginos have masses Mj < 1 TeV. We include transport coefficients for 
bottom and tau Yukawa interactions, recently shown to play an important role in MSSM 



baryogenesis 42]. 



The BMSSM CP-violating sources can have a large impact on baryon number generation, 
illustrated in Fig. Here, we plot ns/s (the baryon-to-entropy-density ratio), normalized 



to the observed value ns/s ~ 9 x 10 11 52j, for maximal AO; i.e., the vertical axis ap 



proximately shows 1/A6 needed to give the observed BAU. In order to highlight the novel 
effects of the BMSSM, we take A(3 = and neglect all MSSM phases (<& = <f> q = 0). The 6 
contributions are suppressed in the large tan/3-limit; we take tan/3 = 3. Other parameters 
are specified in Table HI1 

In the left panel of Fig. |2j we show the BAU induced via 9 contributions to the top (blue 
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dashed) and stop (green dashed) CP-violating sources, as a function the left-handed stop 
mass parameter mg 3 . Since the right-handed stop is light, the stop source is enhanced for 
smaller values of mg 3 due to the resonance of the CP-violating source. The stop source 
becomes suppressed off- resonance, for large mg 3 . The top CP- violating source does not 
depend on rriQ 3 . However, the top-driven contribution to rig/s is suppressed at small values 
of mg 3 due to (i) an enhanced stop contribution to the relaxation rate 4(3], and (ii) more 
charge equilibrating into left-handed stops, rather than tops, reducing the fermionic charge 
available for sphaleron conversion. 

In the right panel of Fig. [2J we show the BAU induced via 9 contributions to the higgsino- 
wino (red dashed) and higgsino-bino (purple dashed) CP-violating sources. Both sources 
are enhanced on-resonance when /i M 12 . In addition, we have studied the bottom and 
sbottom CP- violating sources (not shown); though similar to the top/stop contributions, 
they are suppressed by (mb/mt) 2 and would be viable only for A8 = 0(1). 

The total baryon asymmetry in the BMSSM is a combination of (i) contributions shown 
in Fig. [21 and (ii) standard MSSM contributions induced through A/3 and the MSSM phases. 
Our key point is that the BMSSM contributions can be large, thereby extending the window 
in parameter space that can provide successful EWBG. In the next section, we investigate 
the EDM constraints on this scenario. 



IV. ELECTRIC DIPOLE MOMENTS 



A. EDMs from BMSSM Phases 



EDM searches are sensitive to the same CP-violating phases that generate the BAU, and 
consequently provide powerful constraints on the EWBG mechanism. Currently, the most 
significant EDM bounds are for the neutron and the thalium and mercury atoms 53l-l55|. 
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TABLE II: The parameters used for Fig. \^ in units of Ge V. In addition, we take At 
tan/3 = 3, and all other soft SUSY-breaking masses to be heavier than 1 TeV. 



A b = 0, 
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(a) (b) (c) (d) 

FIG. 3: Examples of BMSSM contributions to CP-violating operators: (a) one-loop and (b) two- 
loop EDM and chromo-EDM, (c) Weinberg operator, and (d) four-fermion operator. 

given at 95% C.L.: 

\d n \ < 3.5 x 1(T 26 ecm , (35a) 

| dn| < 1.1 x 1(T 24 ecm , (35b) 

\d Hg \ < 2.9 x 1(T 29 ecm . (35c) 

In many scenarios, including in our present work, the electron EDM d e provides the dominant 
contribution to g?ti, given by g?ti — — 585 c? e . Under this assumption, the corresponding 
bound is I dJ < 1.9 x l(T 27 e cm [95% C.L.l. 



CP violation in the BMSSM generates, below the weak scale, several classes of CP 



violating, non-renormalizable operators, which in turn give rise to the above EDMs 56j . In 
Fig. EH we show examples of BMSSM contributions to these operators. At dimension five, 
there are EDM and chromo-EDM operators, arising at one- loop order (Fig. [3^l) . Two-loop 
contributions (Fig. [3b) become dominant when first and second generation sfermions are 
heavy (my > 1 TeV) 5 . At dimension six, there are the Weinberg operator (Fig. [3b) 57j and 



four-fermion operators (Fig. [3bl). Novel BMSSM contributions to these operators arise from 
(i) explicit factors of en in the mass terms for sfermions and higgsinos [Eqs. (I17U19j) ]. (ii) tree- 
level scalar-pseudoscalar neutral Higgs mixing [Eq. (1121) ]. and (iii) the complex Higgs vev, 
also arising at tree-level [Eq. (fl6l)]. These contributions give rise to irreducible EDMs that 
cannot be universally suppressed without also destroying the viability of EWBG (barring 
fine-tuned cancellations). 



5 In the BMSSM, one must not push the first and second generation sfermion masses above the scale M, 
where the BMSSM ceases to be valid. In principle, one could consider a modified version of the BMSSM 
in which these sfermions are integrated out along with the physics responsible for the BMSSM operators. 
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We evaluate these EDMs following Ref. 



. This treatment, utilized for the MSSM with 



a loop-induced CP-violating Higgs sector, includes all four classes of contributions shown 
in Fig. [3] and allows for the inclusion of the BMSSM effects described above. However, our 
EDM results are subject to two main theoretical uncertainties. First, the two-loop EDM and 
chromo-EDM contributions are incomplete; the remaining known MSSM contributions 



59 



611 ] have not been generalized to include scalar-pseudoscalar Higgs mixing and cannot be 



easily adapted to the BMSSM. Second, there exist 0(1) uncertainties in the hadronic inputs 
needed for the evaluation of dn g and d n . For the neutron, in particular, there exist three 
different methods for computing d n , each sensitive to a different linear combination of CP- 
violating operators. Below, we show only neutron EDM bounds computed using the more 
recent "QCD Sum Rules" method, which is sensitive to first generation EDM, chromo-EDM, 
and Weinberg operators 62 ]. 



As we show below, the mercury EDM provides the strongest bound on EWBG in the 
BMSSM. (In contrast, EWBG in the MSSM is constrained by d e and d n and is largely 
insensitive to du g .) This scenario will be decisively probed by the combination of future 
EDM searches, which are expected to reach sensitivities of 10 -28 e cm for d n , 10~ 29 e cm for 
d e , and 10 -29 e cm for the EDMs of the deuteron {dn) and proton 63]. 



B. Constraints on BMSSM Baryogenesis 

In this section, we show how limits on EDMs constrain BMSSM baryogenesis. Although 
both MSSM and BMSSM phases, listed in Table fl] can impact both EWBG and EDMs, 
we choose to highlight the BMSSM by setting all MSSM phases to zero. In this case, CP 
violation is governed by the parameters en and e^. 

BMSSM baryogenesis can be driven by squark, quark, or higgsino CP-violating sources, 
discussed in Sec. IIIII In Figs. H] and El we illustrate how current EDM constraints impact 
each of these EWBG scenarios, with parameters given in Table IIIII In each panel, the gray 
bands show the region of the en-e2i parameter space consistent with generating the observed 
baryon asymmetry. BMSSM baryogenesis is inconsistent with en = e2i = (shown by the 
cross), since clearly CP violation is required to generate the BAU. The width of this region 
corresponds to the following range for the wall velocity v w and Higgs vev at the critical 
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FIG. 4: The viable region for stop-driven EWBG in the en-e^i parameter space (in gray), and 
EDM constraints from dn g (red solid), d n (green short-dash), d e (blue long-dash), and do (brown 
solid). Left panel: current EDM constraints; Right panel: future EDM constraints: \d n \ < 10 -27 e 
cm, \d e \ < 5 x 10~ 29 e cm, \do\ < 10~ 28 e cm. (Note: zoomed-in scale on right.) Other relevant 
parameters are specified in Table 17771 



temperature v c (c.f. Sec. IIIII) : 

O.Ol <v w < 0.4, 70 GeV < v c < 110 GeV . (36) 

The edge of the EWBG region closest to the origin corresponds to the minimum value of 
{en — sin2/3e!j) consistent with the BAU, achieved when EWPT parameters fortuitously 
maximize the CP- violating sources (v w ~ 0.03 and such that the v c is largest). 

In Fig. H]we consider the stop-driven EWBG scenario. In the left panel, the regions con- 
sistent with current mercury, neutron and electron EDM constraints (95% C.L.) lie between 
the red (solid), green (short dash) and blue (long dash) curves, respectively. Each EDM is 
consistent with en = €2% — 0. Since the EWBG and EDM bands overlap, this scenario is 
viable in the BMSSM. (In the MSSM, it is not viable, since mg 3 ^> 1 TeV (J.) In the right 
panel, we illustrate how future improvements in EDM sensitivities can exclude this scenario, 
assuming null results. Again, we show the same EWBG band, now zoomed in. The EDM 
bounds correspond to \d e \ < 5xl0 -29 e cm (blue long dash), \d n \ < 10~ 27 e cm (green short 
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FIG. 5: The viable region for top-driven (left) and higgsino-wino driven (right) EWBG in the eu- 
€2i parameter space (in gray), and EDM constraints from dn g (red solid), d n (green short-dash), 
and d e (blue long-dash). Other relevant parameters are specified in Table [7771 

dash), and \du\ < 10 -28 e cm (brown solid). EWBG is excluded within the intersection of 
these bounds. Furthermore, these are pessimistic limits compared to actual expected future 
experimental sensitivities (described above). Of course, the alternative, that EDMs will be 
discovered, is a much more exciting prospect. 

In Fig. we consider EWBG scenarios driven by a top CP-violating source (left) and 
higgsino-wino source (right). As above, viable EWBG can occur in the gray region, while 
the current EDM constraints are shown by the colored bands, as in Fig. [5j 

Strikingly, the bands for EWBG and EDMs in Figs. H] and [5] appear to align in the eu-e 2 i 



EWBG scenario 
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TABLE III: The parameters used for Figs.\J^and\£&in units of GeV. We take tan (3 = 3, e± r = —0.05, 
€2 r = 0.05, Aj = 0, and m 2 ~ = (1 TeV) 2 for all other sfermion masses-squared. 



21 



plane, signaling that they depend largely on the same linear combination of e± t 2i- The BAU 
is proportional to A9, the variation of the Higgs phase across the wall, while the EDMs 
are dominated by contributions proportional to m\ A (the h -A mixing parameter) or 9 
(the vacuum Higgs phase). Indeed, all three parameters are approximately proportional to 
(2eij - sin2/3e 2i ). 

The mercury EDM limit provides the strongest constraint on EWBG in the BMSSM. 
Here, the dominant contribution to d-^g is the down quark chromo-EDM operator induced 
at two-loop (Fig. [3b) with a top quark loop and light Higgs exchange, proportional to 
the h -A mixing parameter m\ A . The subdominant contribution proportional to m 2 HA is 
responsible for the slight skew between the c?H g and EWBG bands in Fig. |H 

Lastly, we describe how EWBG and EDMs depend on the parameters rriA and tan /3. The 
consistency between EDMs and EWBG is generally insensitive to tua- Larger values of tua 
suppress A9 and the baryon asymmetry as ~ l/m A , requiring larger values of ei^j, but the 
EDMs bounds are correspondingly weakened as well. On the other hand, for larger values 
of tan /3, the BMSSM CP- violating sources are either constant or suppressed, while EDMs 
become enhanced — thus leading to stronger constraints on this scenario. 

V. CONCLUSIONS 

Adding dimension-five terms to the Higgs potential of the supersymmetric standard model 
alleviates the fine-tuning problem related to the LEP lower bound on the Higgs mass, and 
has interesting consequences for the electroweak phase transition and dark matter. 

Here, we have investigated the consequences of the new CP-violating phases of these 
terms for supersymmetric baryogenesis and electric dipole moments. Our main observations 
and conclusions are the following: 

1. The introduction of the e\ and e 2 terms implies two new physical CP-violating phases. 

2. Unlike the MSSM, the BMSSM allows for spontaneous baryogenesis, that is baryoge- 
nesis that is generated by a complex phase in the Higgs VEVs that is changing across 
the bubble wall. 

3. In addition, several CP-violating sources that are ineffective in the MSSM, can become 
effective in the BMSSM, namely the stop and top sources. 
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4. It is possible to have successful baryogenesis with all the MSSM phases put to zero, 
and with either or both of the two new phases. 

5. The EDM constraints can be satisfied if the new phases are of order 0.1 or smaller. 
Successful baryogenesis requires, however, that they are not much smaller than the 
EDM upper bound. Thus, barring cancellations, BMSSM baryogenesis predicts that 
EDMs should be discovered if the experimental sensitivity improves by about an order 
of magnitude. 



Appendix A: Higgs phase vacua 

In this appendix we discuss the vacuum value of the Higgs phase 9. At tree level, the 
part of the Higgs potential which depends on 9 is 

\e2\v 2 sin 2/3" 



AVb = -v sin 2/3 Re 



jo 



2|e 1 |^ 1 )- 



(Al) 



At fixed values of v and tan /3, the vacuum structure in the 9 direction is determined by a 
single complex parameter. To see this, define the following quantities: 





= \b\ 


+ 2 






= arj 


l[beS 


], 


|A| 


= \b c s\v 2 


sin 2/3, 


K 




b e ff 


v 2 sin2/3e 4{,?2 


9 


= 9 + As- 



The potential and the minimum equations can be written as 



AV (9) = —\A\ Re 
da AVb (0) = +|A|Im 
d$AV (8) = +\A\Re 
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(A3a) 
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Eq. flA3p admits an analytic solution, depending only on k. The solution is obtained noting 



that Eq. ()A3b|) can be cast as a quartic equation ^2 i=0 Oj sin* 9 = , with coefficients a 



Ihi[k] 2 , ai = 2Im[/«] , 02 = 1 — 4|k| 2 , 03 = — 4Ihi[k] , = 4|k| 2 . Out of the four roots of the 
quartic equation, at most two represent a local minimum of the potential. 
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In order for a double well potential to arise, the e 2 term must be sizable in comparison 
with 6 e g. In particular, for \k\ < 1/2, only one minimum exists regardless of the phase of n . 
The condition on k can be cast as a condition on the mass of the charged Higgs 6 and on e 2 , 

m 2 H — m^r \ ( m 2 H — 

> 3|e 2 | — 2e 2r 5 (condition for single minimum) . (A4) 

2v z 1 \ 2v 2 



The condition (1A4|) depends on neither tan/3 nor e± . For |e 2 j jr | < 0.1, it is fulfilled for 
m H ± ^ 170 GeV. In case that e 2r > , it is enough to impose rriH± > 130 GeV. In fact, 



violation of Eq. (IA4D necessarily implies m 2 H ~ e 2 v 2 |65|. This can be understood in 
general from the Georgi-Pais theorem 66j. (In our case, nonrenormalizable operators should 
be considered instead of quantum corrections to break the CP-symmetry of the potential.) 
Such low values for m# ± (and hence also tua) are, in general, in tension with both direct and 
indirect experimental constraints. We did not pursue further the analysis of this parameter 
regime, even though it may have interesting consequences for EWBG via transitions between 



different phase vacua [67]. 



As a final comment, note that keeping the term of order e 2 in Eq. (IA4D is required for 
small rriH± , where (m 2 H± — m 2 v )/2v 2 becomes a small parameter 0(0.1) . As is the case 
in several points in this work, neglecting independent dimension six operators is still a 
consistent procedure. 



Appendix B: Quantum Corrections 



Quantum corrections to the Higgs sector, in the presence of explicit CP violation in the 



MSSM, were discussed in detail in Refs. 
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681 ] for zero temperature and Refs. 



29 



30j for 



finite temperature. Here we present only the corrections that are directly relevant for the 
vev phase in the BMSSM. We consider squark, chargino, neutralino and scalar Higgs loops. 

We write the potential as V = Vo+AVi+AVr, where V is the tree-level, zero-temperature 
part, given in Eq. (j4j), AVi is the zero-temperature one-loop part, and AVr is the finite- 
temperature correction. For AVi, one sums over all particle species with field dependent 



6 For k on the real axis, the emergence of an additional solution with 6^0 implies spontaneous CP 
violation. This situation was discussed in 



64| for the BMSSM. 
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masses: 



AV 1 = J2 nim - ) 



64vr 2 V'" Q 2 2; ' (B1) 
where Q is the renormalization scale, which we choose as Q = m t . For AVt, we include the 
pressure and daisy terms: 

'm 2 (0)\ \ - rtiT 
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where the second sum includes scalars and longitudinal gauge bosons and with the J func- 
tions defined by 



dyy 2 \n(lTe-^ 
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Beginning with zero temperature, it is useful to first obtain an analytical estimate of the 
contributions from stops, charginos and neutralinos. We follow the procedure introduced in 
Section [Til an d apply it to Vo + AVi. We find the corrected expressions for \b\ and 9, 
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The quantity 5V(oc dV/d6) encodes the various contributions: 
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where, to 0(ei) and 0(g 2 ), and neglecting contributions to the charged Higgs mass arising 
from diagrams involving chargino-neutralino and stop-sbottom loops, we have 
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The loop function is 
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The conditions that the CP violation induced by the usual MSSM loop corrections be- 
comes comparable to the tree level nonrenormalizable contribution can be written as 

ei 3M 2 /i ei M 2 fJ, 



(7 2 /167r 2 ~ 2v* 0.05 ~ (600 GeV) 2 ' 

e x < 3A t fi ex A t \i 



y2/167r 2 ~ 2t> 2 0.05 ~ (400 GeV) 2 ' 

We learn that the BMSSM terms can easily dominate. Furthermore, if <px,2,t <C 1, the 
radiative corrections are unimportant as they mostly serve to slightly shift the value of |6|. 

Eq. (IB6P includes also a contribution from ex itself, arising via its appearance in the squark 
and sfermion mass matrices. Here, the stop contribution is negligible compared to that of 
charginos and neutralinos. This reflects the fact that ex corrects only the mass splitting of 
stops, while it enters the trace of higgsino mass matrices. The ex loop correction tends to 
cancel the tree level term. It becomes relevant if fi 2 /v 2 is large enough to compensate for 
the loop suppression; in practice, the term is significant for jj, > 500 GeV. 

Finally, a common feature of the zero-T corrections of Eq. (IB6I) is that they depend only 
mildly on the vev. Hence, the net effect of these terms in the vicinity of the origin of field 
space is to globally shift the value of 9. This implies that the phase variation across the 
bubble wall, AO, and consequently the novel BMSSM contributions to the BAU, remain 
mostly unchanged. 

Proceeding to finite temperature, let us again obtain some analytical understanding, 
beginning with the stop sector. We consider the plausible limit where the heavier stop, t 2 , is 
Boltzmann suppressed, while the contribution of the lighter tx admits a high-T expansion. 
Then, 

M0,T) ~ |6|+2|e 1 |0V ?1 + ^ r (|A t/ i|e*+2|e 1 |0 2 c 2 e^) 

,2rp2 



Sy^Atfiex 

While there is no loop suppression in Eq. flB8j) . the thermal correction at the critical temper- 
ature is still down by a factor ~ T c 2 /mg 3 . Note also that the MSSM term oc |^lt/x| e 2<?it is field 
independent to O (T c 2 2 /mQ 3 ) , such that its contribution to the variation of 9 across the 
bubble wall is suppressed. Using Eq. (1B8|) . we can compare the tree level effect of the non- 
renormalizable BMSSM terms with the leading thermal correction of the MSSM. Focusing 
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on the phase variation along the wall, 
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In the bulk of this paper, we isolate the novel BMSSM effects by using small or vanishing 
values of A t , as well as vanishing <fi t . Hence only the term proportional to ei0 2 c| in Eq. ( TB8j) 
remains. Since this term is field dependent, it does not affect the determination of 9 in the 
symmetric outskirts of the bubble wall, where <p <C T c . Since it is doubly-tan /5 suppressed, it 
can typically be neglected with regard to the variation of 9 along the wall, even for niQ, A ~ T c . 

Moving on to the Higgs and higgsinos, we find that Higgs-Higgs and Higgs-higgsino 
interactions lead to a non- negligible shift in the finite-temperature Higgs phase 9{z) at 
0(ei). This effect is easy to understand by considering a high temperature regime T > m, 
where m denotes any mass parameter in the Higgs and (weak) gaugino sectors. (Of course, 
at such high T > T c the universe is globally symmetric. Nevertheless, it is illuminating to 
provisionally pursue this line of argument.) In this temperature regime and near the origin 
of field space, the one-loop thermal potential due to Higgs and higgsino particles can be 
expanded, 

rp2 J,2 rp2 rp2 

Vt' R ~ - TrM^ o + - TrA4 ± + - TrM^ + - TrXtX 

D - 3|£l|T 2 202g2/3 cos(e + A), (BIO) 



The phase- dependent piece of the potential becomes 
V T {6) ~ -0 2 s 2/ 3 



\b\ cos 9 + 2|ei| ( 2 + — \ cos(6 + ^) - cos ( 2 ^ + # 2 ) 



.(Bll) 

4 J 2 
In the high-T approximation and in the cot /3 ^> \ei\v 2 /m? A regime, setting T = T c , the 
values of 9 in the symmetric (s) and broken (b) phases are 

9 S = hm 9{z) ~ -arctan ' l ± c B12 

2->-oo 2\b\ + 3 1 e x IT/ cos 17 1 rn z H± s 2 ^ 

9 b = lim 9(z) ~ 6 S + A6 . (B13) 

z— >+oo 

In other words, the dominant effect is to shift 9{z) by 9 S over the tree-level result [Eq. f[2T|) ]; 
however, the relative phase (#& — 9 S ) that governs EWBG remains unchanged from its tree- 
level value of Ad. 

We should make the following comment, regarding the calculation of 9 S in Eq. (1B12I) . 
Any dependence of the potential on 9 must be proportional to (f) 2 s 2 /3. Thus, in the truly 
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symmetric region where the vev vanishes, the value of 9 is not well defined. The quantity 
we denote by 9 S corresponds to the value of the CP- violating phase in the almost symmetric 
regime, where the vev is finite, but smaller than any other mass scale in the problem. The 
non-zero value of 9 S indicates that, in the complex plane, the origin can be approached from 
different directions. Since 9 cannot affect the dynamics in the symmetric regime, we may 
extend the definition of 9 S to z — > — oo. Having clarified this point, it is important to keep in 
mind that in the region of <p < gT there are additional non-perturbative corrections to the 
scalar potential, which we do not consider here beyond the daisy resummation introduced 
in Eq. (IB2I) . We note that the daisy corrections, which include ^-dependent terms in the 
scalar self-energies, do not significantly affect our results for 9 S . 

In Fig. [HI we illustrate the role of quantum and thermal effects in dictating the complex 
vev phase at zero and finite temperature. We go beyond the approximation of Eqs. flB12j) 
and (1B13P by using the full one-loop thermal potential, instead of the high-T expansion 
given in Eq. ( 1B10I) . We solve for the asymptotic values of the complex phase by minimizing 
the potential at the critical vev and near the origin of field space, setting T = T c . Since 
the ei loop contribution depends mainly on /i, we use fi as an independent variable. To 
emphasize the role of the T-dependent terms, we fix all parameters and repeat the plot in 
two panels, once for T c = v c = 100 and once for T c = 150 , v c = 110 GeV. We minimize the 
potential numerically, accounting also for the zero- and finite-temperature effects of the Higgs 
scalars. Two main features which were mentioned above are worth pointing out. First, the 
onset of the zero temperature effect appears at large \x ~ 400 GeV. Second, while quantum 
corrections shift 9 S sizably, the value of A9, the phase difference between the broken and 
symmetric domains, is much less affected. 
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FIG. 6: The vev phase, including quantum and thermal corrections. Unspecified parameters are as 
in Table\U Left: T c = 100 GeV, v c = 100 GeV. Right: T c = 150 GeV, v c = 110 GeV. Green, red 
and blue curves correspond to 9 S , 9^ = 9 S + A9 and the zero-T value of 9, respectively. Dashed 
(solid) curves denote tree level (one loop) results. 
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